CUDLEV  .SCHOOL 


AVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


Implementation  of  a  Reliability  Shorthand 
on  the  TI -59  Handheld  Calculator 

by 

Hans -Eberhard  Peters 

October  1982 

Thesis  Advisor:  J.D.Esary 


Approved  for  public  release;  distribution  unlimited. 


1JT205425 


SECURITY   CLASSIFICATION  OF   THIS  PAGE  (Whan  Dm  Entatoa-) 


REPORT  DOCUMENTATION  PAGE 

1   REPORT  NUMIIR 


2.  OOVT  ACCESSION  NO 


4.   Ti  tl£  and  Subtttl*) 


Implementation  of  a  Reliability  Short 
hand  on  the  T I  -59  Handheld  Calculator 


7.     AuTwOR.«> 


Hans-Eberhard    Peters 


9.     PERFORMING  ORGANIZATION   NAME   ANO   AOORESS 

Naval  Postgraduate  School 
Monterey,  California  9  3940 


I  I.     CONTROLLING  OFFICE  NAME    ANO   AOORESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


TT     MONITORING   AGENCY  NAME  *    AOORESV"  dlllatant  tram  Controlling  OlUea) 


Naval  Postgraduate  School 
Monterey,  California  9  39  40 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3       RECIPIENT'S  CAT  ALOG   NUMSEl 


S       TYPE   OF    REPORT    4    PERIOD   COVERED 

Master'     Thesis ; 
:tober    1982 


9^ 

•■     PERF 


ORMING  ORG.    REPORT   NUMBER 


•  •     CONTRACT   OR   GRANT 


NT    HUltRll, 


10.  program  Element,  project  task 
AREA  a   WORK  UNIT  NUMBERS 


12.  REPORT  OATE 

October   1982 


II.     NUMBER  OF   PAGES 

72 


IS.     SECURITY   CLASS,   (ol  thf  roport) 

Unclassified 


'»«.     DECLASSIFICATION/  OOWNGRAOlNG 
SCHEDULE 


16.     DISTRIBUTION   STATEMENT  :ol  Ihlo  Mapottj 


Approved    for   public    release;    distribution    unlimited. 


17.     DISTRIBUTION  STATEMENT  (ol  tho  aaatracl  antarad  In  Block  20.   II  dltloront  tram  Roport) 


IS      SUPPLEMENTARY  NOTES 


19.     KEY  WORDS  (Conllnuo  on  roworoo  tlda  II  nacoooary  and  Identity  my  aloe*  ntombor) 

TI-59 

Handheld    Calculator 
Programmable    Calculator 
Reliability    Shorthand 


20.     ABSTRACT  (Continue  an  ravaroa  old*  II  naeaaoawy  and  Idanllty  ay  oioeM  numaar) 

It  is  shown  how  a  reliability  shorthand  can  be  implemented 
on  a  handheld  calculator. 

Assuming  constant  failure  rates,  basic  structures  are 
used  to  show  how  the  shorthand  can  be  applied.  Several 
examples  are  worked  out  that  show,  how,  with  component 
failure  rates  as  input,  a  handheld  calculator  can  be 


DO 


FORM 
I    JAN    7 J 


1473  EDITION   OF    1   NOV  •»  IS  OBSOLETE 

S/N    0  102-014-  660  1 


IECURITY  CLASSIFICATION  OF  TNIS  PAGE  (Whan  Data  Sntarad) 


^cuw*  Ci  »«ti>'C*TtOM  py  Twit  mi^w  n...  !-<•«•« 


DD      Form       1473 

1  Jan  73 
S/N    0102-014-6601 


tccu«t*v  eL*Mi*t<ArioN  0'  T*»«  **a*f*»«"  a«««  f«#»»»*i 


Approved   for   public   release,    distribution   unlimited, 


Implementation   of    a   Reliability   Shorthand 
on   the   TI-59    Handheld    Calculator 


by 


Hans-Eb erhard   Peters 
Major, German    Air    Force 
Dipl.-Betriebsv. ,    Fachh cchschule    des    Heeres    I,    1974 

Submitted   in  partial   fulfillment   of   the 
requirements    for    the   degree   of 


MASTER    OF    SCIENCE    IN    OPERATIONS    RESEARCH 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
October,  1982 


"ONTER^y^^^-SCHOOL 


ABSTRACT 

It   is   shown    how   a    reliability    shorthand   can    be   imple- 
mented on   a   handheld   calculator. 

Assuming   constant    failure  rates,    basic   structures   are 
used    to   show   how   the    shorthand   can    bs    applied.    Several    exam- 
ples   are    worked    out   that   show,    how,    with    component    failure 
rates   as   input,    a    handheld    calculator    can    be   used    to    compute 
the    reliability   of   a    sys-em. 

Two    TI-59    programs   are    provided    2.3    a    computational   aid. 
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I-   INTRO  DUCT  ID  N 

Systems  and  components  can  be  in  either  of  two  states: 
either  they  are  functioning  or  they  have  failed.  The  abil- 
ity, that  a  system  stays  functioning  over  a  predetermined 
time  interval  is  called  its  reliability.  It  is  generally  not 
realistic  to  assume  that  a  system,  say  a  lightbulb,  will 
fail  at  a  specified  time,  but  rather  that  T,  the  time  to 
failure,  is  a  random  variable  which  has  a  probability  dis- 
tribution that  can  be  specified.  The  probability  distribu- 
tion for  a  time  to  failure  is  called  its  life  distribution. 
In  this  paper  we  will  solely  be  concerned  with  one  specific 
type  of  life  distribution  which  is  especially  important  in 
reliability  theory  and  practice,  the  exponential  distribu- 
tion. It  has  the  property  that  the  remaining  life  of  a  used 
component  is  independent  of  its  age  (the  "memcryless"  prop- 
erty) ,  i.e.  a  functioning  component  is  always  as  good  as 
new,  the  failure  rate  is  constant.  The  memoryless  property 
is  the  basis  for  a  reliability  shorthand,  one  that  can  be 
implemented  on  a  handheld  calculator. 

Depending  on  the  size,  structure  and  life  distribution 
of  a  system,  probability  statements  about  its  time  to 


failure   are   in   general   not    easily   achieved.    Forming   the   sum 
of    independent    life    lengths     (i.e.      convolving   tne    corres- 
ponding  life   distributions)     requires    knowledge   of    integral 
calculus   and  computations  can  become    rather   tedious. 

In   the    case   of   the  exponential   distribution,    though, 
computations    can    be    simplified   by   translating   the    problem 
into    a  simple    shorthand   notation   and    using    this    shorthand   as 
input    for   some   computing   device. 

In   this   paper    we    will   show  how   a    reliability    shorthand 
can    be   implemented   on   a   handheld   calculator.    Basic   struc- 
tures   are   used   to   show  how    the   shorthand   can    be    applied.    Two 
TI-59    programs    are    provided    as    a   computational   aid.    Formulas 
for    the   convolution   of  up  to   four   exponential   random    varia- 
bles   can   be   found    in   Appendix  A.    Appendix   B   contains    a   user 
guide    to   the   TI-59    programs. 


II.       THE    CONCEPT    3F    A    RELIABILITY    SHORTHAND 

A.  BASIC    NOTATION 

The    survival   function    of  a    life    length    can    be    derived 
from   the   distribution    function. 
Let 

T  :    life    length 

F(t)  =  P(T^t)  be  the  distribution  function  of  T 
Then 

F  (t)  ■  P  (T>t) 
=  1"F(t) 
is  the  survival  function  of  T. 

In  the  case  of  the  exponential  distribution,  F(t)=  e    , 
where  A.  is  the  failure  rate.  Translated  into  shorthand,  the 
life  distribution  is  denoted 
EX?  (  I    )  . 

B.  CONVOLUTION  OF  DI STRIBUTIONS 

When  independent  random  lives  are  summed  up,  the  corres- 
ponding life  distributions  have  to  be  convolved  to  determine 
the  probability  that  the  sum  of  the  lives  will  exceed  a  spe- 
cified  time  t.      Let 

T,  ,  T2  :    independent    life    lengths 

10 


'hen 


F,  (t)  , F2  (t)   :  the  corresponding  survival 

fa  notions 
f*  (t)  ,fi  (t)   :  the  corresponding  density  functions 
T  =  T^  +   Hz      :  the  total  life  length 

F(t)  =  P(T>t) 


=  P(T^  +  T2  >t) 

t 

=     %     (t)   +  j?2    (t-S)  f,  (3)  dS. 
0 

This  means  that  T  will  exoeed  a  specified  time  t  when 

-either  1A    exceeds  t 

-or  T.,  is  smaller  than  t,  say  equal  to  s,  and  T4 
exceeds  t-s. 
Integration  with  respect  to  s  (i.e.  summing  over  all  possi- 
ble values  of  s  )  is  called  the  convolution  of  T.,  and  Tz . 
When  TA    and  T2  are  both  exponentially  distributed  with  fai- 
lure rates  K^     and  /t2  ,  i.e. 


F„  (t)  = 


-lot 


F2  (t)  =  e 


-lx-t 


then    the   survival   function    of  I    is 


-M      ?-iz(*'*)j   :^s 


F  (t)    =    e  '***    +   Je" 


ds. 


Translated  into  shorthand,  the  survival  function  is  denoted 
EXP  {i.)     +  EXP  (  I,   )   . 
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This    shorthand    notation    is    heuristically    apparent.    We   can 
visualize   a    1    component   /   1    spare   system    with    EXp(/t4)    and 
Exp  ( /i2  )    lives   respectively.      From   component    1    the   system 
has    an    EXP  ( LA  )     life    to   begin    with.    When    component    1    fails, 
the    system    has   an   extra   EXP  (  <lz)    life. 


MIXTURE    OF    DISTRIBUTIONS 


1.      MIX-Notation 


In    the    previous   chapter,    we    formed   the   sum    of    inde- 
pendent  random   lives,    which    each    had    weight    one,    i.e. 


T   =   T,,    +    T2  . 


Now    consider 


r 


T,,      with   probability    p. 


T   =     < 


Tz      with    probability    p£ 
where    p.,    ♦    p2   «    1    . 

Let    D^    and    D2  be   the    probability   distributions   of    the   random 
variables   T1  and  Tj,    respectively.    The    corresponding   survival 
functions   are   F,,  (t)    and   F2  ( t)  . 
Than 

F(t)     =    P„F~    (t)  +    p,  F,  (t)  . 


i    1 


2    -2 


In  shorthand,  the  mixture  of  distributions  D^  and  Dz    with 
respect  to  the  mixing  probabilities  p„  and  p,  is  denoted 
MIX  [  P,  D^   ,  p202  ]. 
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2.      Distributive    Law 


Now   let 


where 


T    =    T3    +    T» 


T,,    with    probability    p 


T»    = 


Then 


Tt   with    probability    1-p. 


r 


T ,,   with   probability    p 


T      =    T3    ♦  < 


Ti    with   probability    1-p. 


r 


T,     + 


3    *    l1 


with    probability   p 


T      =    < 


Tj    +  T2      with    probability    1-p. 

The    distributive    law    holds    due   to    the    fact    that    the    sura   of 

the    mixing    probabilities    for   T>,    and    r2    is    equal    to   one. 

The    survival   function    of    T    can   be    found    by   convolution: 

•fc 

F(t)    =    F5  (t)    +    J(pF4(t-s)     ♦     (1-p)F2  {t-s))f3  (s)ds. 

o 

With  D,,  ,  D2  ,  Dj  being  the  probability  distributions  for  Tn, 

T2  ,  Tj ,    the  distributive  law  can  be  applied  to  the  shorthand 
notation: 

D3+  MIX  [pD,,,(1-p)D2]  =  MIX  [  pUD,,  +  D3  )  , 

(1-p)  (D2  +  D3)  ]. 

13 


Graphically  this  can  be  represented  as  follows: 


D3  + 


D, 


D. 


1-D 


D-,  +  D3 


D2  +  D3 


Figure    1:      Distributive    Property   of   the    BIX- Notation 

3 .      Degeneracy   at    the   Origin 
Let 

P(T=0)     =    1. 
Then    the   distribution    of    I    is  degenerate    at    zero. 
In    shorthand   notation,    such    a   distribution    is    called    the 
ZERO- distribution. 
Now    let      T    =    T,,    +    10 

where   T4    and   T0    have    probability    distributions    D,,     and    ZERO 
and    survival   functions  F^  (t)    and   f0  (t)    respectively. 

Then 

t 

F(t)  =  f;  (t)  +  jf0  (t-s)f.,  (s)ds 
=  F«  (t ) . 
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The    ZERO-distribution    doesn't   add   anything   to   another   dis- 
tribution,   so   for   instance 

D^    +    ZERO    =    D^ 

D2    +    MIX[pD^  ,(1-p)ZERO]    =    MIX[p(D,,+    D2)  ,  (1-p)  Dg  ], 
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III.       APPLYING    A    RELIABILITY    SHORTHAND 
After   this    brief    survey    over   the   concept  of   a   reliabil- 
ity   shorthand   we   will    now   show   how   the    shorthand    can    be 
applied.    To   do   sc    we    will   use   basic    structures.    Part    A   of 
this    chapter    will    give  examples    whose    representation    in 
shorthand   requires   only    basic  notation    described    in   Chapter 
II,    Parts    A  and    B,    whereas    Part    B   of   this    chapter    will    give 
examples   whose   representation   in   shorthand    makes    use    of   the 
MIX-notation   and  the    ZERO-i istr ibution. 

A.       SUMS    OF    EXPONENTIALS    WITH   WEIGHT    ONE 
1 •      Simple   Series   Sy_st e m 

A   series   system   is    a   system    which    is    functioning, 
when    all   its    components   are    functioning.    A    two-component 
series  system   can    be    graphically    represented   as    shown    in 
Fig. 2    . 
Let 

T    :    life    of    the   system 
Ti  :    life    cf    component    1 


T^  :    life    cf   component    2 

F4  (t)     =    survival    function    of    component    1 

— -        r\ 
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Figure    2:      Two-Component   Series   System 

F2  (t)     =    survival    function    of    component    2 
=   e      i 
Then 

T   =   min  (    T4  ,    T2  ) 

F(t)    =    survival   function    of    the    system 
=    P(    min     (T4  ,    T2)     >   t) 
=    P(    T<>    tf     T2>    t) 
Assuming   independence    of    the    two   components 
F(t)    =    ?(    T^>    t)     P(  T2>   t) 
=    R,   (t)     F2  (t) 
=    e  e 


-U, +**)  £ 


=   e 


The    shorthand   notation   for    this    sys-sm    is 
EXP     (/{„  +    42)  . 
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This    is    intuitively   apparent,    as   the    system    has   an    exponen- 
tial  survival   function  with    failure   rate  A*A  +    Lz 
2  .      Si nrple   Parallel    Sy_s  tern 

A   parallel   system   is   a   system    which    is    functioning, 
when    at    least    one    of    its   components   is    functioning.    A   two- 
component   parallel   system   can   be    graphically   represented   as 
follows: 


Figure    3:       Two-Component    Parallel    System 


Let 


Then 


T,,  ~  EXP(   i    )     ,    T1  ~EXP(  I    ) 

T   =    max  (T,,,    Tz) 

F(t)    =    P(max(T,,,     Tt)    >    t| 

=    1-    P  (max(T^  ,    T,)    ^   t) 
=    1    -    P(I^    tf    T2  <    t) 
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Assuming  independence  of  the  two  components, 
F(t)  =  1  -  P(T^  t)  ?  (T2£  t) 

=    1    -    *,  (t)      F2  (t) 

-At  -vl* 

=    1    -    (1-e*T)  (1-e        ) 

„       -*t  -2A.t  t 

=    1    -    (1-2e  +    e  ) 

~4*  -2  It 

-   2e  -   e 

The    shorthand   notation   for    the   system    is 

EXP  (2  >i  )     +    EXP  (  K  )  • 
This    follows    intuition  as   the  system    has    an    EXP(2/i)     life    to 
begin    with    and   when   one   component   fails    it    has  an    extra 
EXP  (  K  )    life   due   to   the    memoryless   property   of   the   exponen- 
tial  distribution. 

3  •      St andtyy^yste  m   H-L%h    Dissimilar    Components 

Suppose   a    system   consists   of    two    components,    one 
active   and    one   spare.      The    active   component    stays    in    service 
until    it    fails    and   then    immediately    is    replaced    by   the 
spare. 

Let    the    time    to    failure    of    the   two    components    be 
T,,  ~  EXP  (/{,,,  )    and   T2  **•  EXP  (  /t2  )    respectively. 
Then    the   system   time    to    failure    is 

T   =   T,,    +    T2 
and   the    survival   function   of   the    system    is 
F  (t)     =    P  (T    >   u) 
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Figure   4:       Standby    System 


=   p 


I 


-a  (t)    +      F2  (t-s)f^  (s)ds 


v 

J- 


e  +    I  e  a,,   e  ds 


•4,  5 


/L      -4»* 


/t„  -  ^j 


/L-2  ~/t<1  * 


A>A  -  A  i 


The    shorthand   notation   for    the   system's    survival    function 
should    be   obvious.      The    system   has    an    EXP(  hA  )    life    from    the 
active   component    and    an    additional   EXP(/tj)    life    from    the 
spare.    So    the    shorthand   notation    is 
EXP  (iA)     *    EXP  (  Az  )  . 
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B.   SUMS  OF  EXPONENTIALS  WITH  WEIGHT  BETWEEN  ZERO  AND  ONE 

The  examples  given  in  the  previous  chapter  only  involved 
exponential  lives  with  weight  one.  Ndw  we  will  look  at  some 
structures,  whose  survival  function  has  a  shorthand  notation 
which  includes  the  MIX-notation  and/or  the 
ZERO- distribution. 

1  •   P^£^Ii£l  System  with  Dissimil  ar  Failure  Rates 

The  notion  of  a  parallel  system  has  been  introduced 
in  Chapter  III. A. 2  .  We  now  look  at  the  case  where 

T,,  ^  EXP  (^  )  and  T2  ^  EX?  (  il  )   . 
Then 

T  =  max  (T,,  ,  T2 ) 

F(t)  =  P(  max(T,f  T2)  >  t) 

=    1-P(    max(T^  ,    Ta)    £    t) 
=    1-P(T,<    tf    T2^    t) 
Assuming   independence    of    the   two    components 
F(t)     =     1-P(T,*    t)      P(T^    t) 
=    1-Fi   (t)     Pj,    (t) 
=    1-d-e"^*  )  (1-6^**) 


=    1-  (1-e  -e  +  e 

=    e  +e  -e 


zj* 
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To  find  the  shorthand  notation  of  the  system  consider  all 
the  ways  which  lead  to  the  survival  of  the  system: 
-either  both  components  survive 
-or  component  1  fails  and  component  2  survives 
-or  component  2  fails  and  component  1  survives. 


If  one  component  fails  and  one  survives,  in 


Jit 


l«+  A: 


fraction 


of  the  cases  the  survivor  will  be  component  1  and  in 
fraction  of  the  cases  it  will  be  component  2. 
This  can  graphically  be  represented  as 

EXP(/l„  ) 


EXPfX,,  +  iz)     + 


I. 


4- +  4: 


iA  +  i 


EXP(  A,   ) 


Making   use    of   the    MIX-notation   the   shorthand   notation    then 


EXP(/t„   +  A,z)    *    MIX[         L      SXP(  An  )  ,T^L-3XP(/i2  )   ] 

and    using   the   distributive    property   it    becomes 


MIX 


[       lli     (EXP(^)     +    EX?(  ^  +   /i2)  )     , 


>t«»     +    /Vj 


/U     +     /t' 


-(EXP  (  X2)     +     EXP  (/{,„+  lz)  )      ]. 
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As    a   check    to    see   that   this   shorthand   notation 
represents    the   survival    function    of   the    system,    we   derive 
the   survival  function    from    the   shorthand   notation: 

— —  (  t  +    J  £  *«£  ^ s  J 


F(t)     =   j 


t 

(«2        +    J£  <i8£         4sJ 


k+  +  At 


e  +9  -   e 


This   verifies    that    the    shorthand   notation    indeed 
represents    the   system's   survival   function. 
2«      Series   System    with    One   S^are 

Let   us   now    look   at    a  two-component   series    system, 
whose    components   have    dissimilar    failure    rates   with   one   com- 
ponent  having    a   spare: 

Component    1    has   ths    constant    failure   rate  X,,    and 
component    2   and    the   spare   have   the   constant    failure    rate  /li  . 
The    spare   can   only    rsplacs    component:    2. 
Let 

F^  (t)     :    the    survival   function    of    component    1 
?z  (t)     :    the    survival   function    of    the   standby 
system   component    2    with    its    spare. 
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Figure   5:      Series    System   with    one   Spare 

The    survival   function    for   a    standby   system   was   derived   in 
Chapter    II.  B.      Therefore 


F,  (t)    =    e 


v  t 


lz 


ds 


=   e"^*     +      i^1^    j  is 

=   (1+  izt)   r^*    . 


F,  (t)    =    e 


-/i.t 


Now 

Then 

(1     +    *2  t) 

To   translate   the   survival   function   into    shorthand    notation, 


F(t)   =  F„  (t)  F2(t) 
=  (1  +  4,  t»  9 


let  us  consider  the  ways  in  which  ths  system  can  survive: 
-either  bcth  components  survive 
-or  component  2  fails  and  its  spare  survives. 
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K. 


-fraction  of  the  time  it  will 


If  one  component  fails,  in 

be  component  1 ,  which  means  that  the  system  will  not  sur- 


vive;   in 


X; 


fraction   of   the   time   the    failing   component 


will    be   component   2. 

This    can   graphically    be   represented    as 


EXP  (  L     +   &i    ) 


ZERO 


lx 


EXP(  A,*  Az) 

Using   the    MIX-notation   the    survival    function    then    is 

EXP  {  AA     ♦    iz    )     +    HIX[- — 7-  ZERO,-; — —-EXP  (  bA  +  iz  )      ] 


=    MIX[ 


L 


=    MIX[ 


I, 

i  +  ^z 

/U 

K 

+  k% 

I, 

k 

1  +  A>z 

il 

(ZERO    ♦  EXP(    I*    +   A*z    ))     , 

(EXP  (  /t,  +  Al  )     +    EX?  (  £4  +  Xi  )     )      ] 

<EXP(  i,  *  lt))     , 

1T-(EXP(  h  +  li  )     +EXP(^+/[2)     )      ]. 


To   prove,    that    the    shorthand   notation    does    represent    the 
survival   function,    we    derive  the    latter    from    the    shorthand: 


4A  ~(A«*A.i)t 

f  (t)  =  - — —  e  + 


/t 


-  (l«-*A,i)-k 
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=   (  1  +  A2t  ) 


-(<L«  +  Az)-t 


This  is  the  previously  found  result  and  this  verifies,  that 
the  shorthand  notation  does  represent  the  system's  survival 
function. 

3  -      Two -out -of- Three    S  y  st  em 

As    a    last    example    in   this    chapter,    we   will    look   at   a 
Two-out-of-   Three   system. 

Consider   a    three   component    system,    whose    components 
have   constant    failure    rates  AA  ,  /i^and   A$   respectively.    The 
system   is    functioning,    as   long    as   two    out    of    three   compo- 
nents  are    functioning    (see    Fig.    6). 

In    other  words,    the    system    is    functioning    as    long    as 
there    is   a    path    through   the    system    . 

Alternatively,    the    system    can    be    visualized   as    a 
par allel -series   system    (compare    Fig.    7)  . 
The    survival   function    of    -he   system    is 

F(t)     =    ?(T^     >t   A    T^    >t)     +    ?  (T,,     >t  A    T5    >t) 
+    ? (T2    >t  A    Tj    >t) 

-  P((T,,    >t  A    ?z    >t)    a     (T,,    >-  A    Ts     >t)  ) 

-  P((T^    >t  A    Tz    >t)    A     (T2    >t  A    ?3    >t)) 

-  P(  (T^    >t  A    P3     >t)    A     (T2    >t  A    T3    >t) ) 
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±A 

■  -          "            "     '            T 

1 .. 

A2 

A3 

-  » 

Figure   6:      Two-out-of-Three   System 


..,,..,.             .               , 

A 

A2 

'     A„ 

A3 

_  , ......  i 

A2 

A3 

Figure   7:      Two-out-of-Three    System 

♦    ?(  (T^    >t  A    T2    >t)   A     <T„    >t  a   T,    >t) 
/\  P(T2    >t  /\    T3    >t)  )  . 
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Thus 

F(t)     =    P(T^    >t  A    T2    >t)     f    ?(Ti     >tA    T3    >t) 
+    P(T2    >t    A     T3    >t) 

-  3P  (T^    >t  A     Tj     >t  A     T3     >t) 
+    P(T/,  >t  A    T2   >t  A    r3     >t) 

Therefore,    and   assuming   independence    of    rhe   components, 
F(t)     =    P(T^    >t)     P  (?-,    >t)     +    P(T,,    >T)     P(T3    >t) 
+    P(T2    >t)     P(T3    >Z) 

-  3P  (T^    >t)     P(T2   >t)     P  (Ts    >t) 
+    P(T^    >t)     P(TZ    >t)     P(T3    >t) 

=    P(T^    >t)     P  (T2    >t)     ♦    P(T^     >t)     P(T3    >t) 
♦    ?(TZ    >t)     P  (T3     >t) 

-  2P  (T4    >t)     P(T2    >t>     P  (T3    >t)      " 

-U„+4»H  -ll«4j.3n  -di+<Lj)t 

=   e  +e  +e 


-    2e 


-(  4<7  *  ki  +  At)  ir 


Now    let    us    consider   all   the    possible    ways,    in   which   the    sys- 
tem   can    survive: 

-  either   all   components   survive 

-  or    component    1    fails    and    component    2   and    3 
survive 

-  or   component    2    fails    and    component    1    and    3 
survive 
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-  or  component  3  fails  and  component  1  and  2 


survive. 


M 


If  a    component   fails   and   the   other   two    survive,    in  — 
fraction   of   the   time    it    will   be    component    i    ,    i    =    1,2,3    . 
This   can   graphically   be   represented    as 


EXP(<1*  +  i%  +  >lj  )     + 


EXP(  /(,,+    X»    ) 


EXP(  iA  +    1}    ) 


EXP  (/{,.,  +-     iz   ) 


The    shorthand   notation   then    is 


EXP  (  lA+  \t  +  &3)      +    MIX[ 


L 


•EXP(  ^  +  /i3)  , 


-EXP(  X„  +  X,) , 


=    MIX( 


X. 


L 

+  /i? 

AiA 

*  <lz 

A., 

+  /i3 

J —EXP  (  iA  *  /i^)   ], 


(EXP(/i2*^,)     +    EXP(/l„  +  /L  +/U  ))  , 


2  T  ^3 


/i/)  +   /tz  ^Aj 


(EXP(^  +  /j)     +    EXP  (  l„  +  lz  +X3  )  )  , 
(EXP(^/|+^2)     ♦    EXP(X„  f  ^  +  /i3  )  )      ], 
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Again,    as    a   check   that   the    shorthand    notation    represents   the 
survival   function,    let   us    derive   the    survival   function    from 
the    shorthand   notation: 


F(t)     = 


I 


A,  A  +    K  X   +  ^  "J 


tc 


-(A,t+A,*H 


1    £ 


X 


/Ia+  Ai  +  A>i 


[«" 


-    llA+lz)t 


-It 


-LIa  +A,z+  A,z)-b 


The   result    again    proves    that    the    shorthand    notation 
indeed   represents    the    survival   function    of    the   system. 
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IV.       IMPLE MEN TING    THE    SHORTHAND    ON    THE    ri-59 
The    concept  of   a  reliability   shorthand    is   introduced    in   the 
course    "Reliability   and   Weapons    System    Effectiveness 
Measurements",    OA    4302,    at    the   Naval    Postgraduate    School, 
Monterey.    Most   students   taking   the   course   are    in    the 
Operations    Research    (OR)     -    Curriculum. 

The    choice    of   the    TI-59    as   the   computing   device,    on 
which    the    shorthand   was   to    be   implemented,    was   based    on   the 
fact,    that    each    student   in    the   OR-Curriculum   is    issued   a 
TI-59    for    use    in    basic   probability   and    statistics    courses. 
Thus    almost   every    student   at   the    Naval    Postgraduate    School, 
who    is   introduced    to    the    shorthand,    is    familiar    with    the 
TI-59   and. has    access    to   such   a   calculator. 

A    program,    that    uses   the    shorthand    notation,    times    to 
failure   and   failure   rates   as   input,    should 

-  calculate   the   survival    probability   of    basic 
structures  /   small   systems   and 

-  require    moderate  computation    time. 

To    achieve    these   requirements   it    was    decided   to    incorporate 
all    solutions    for   the    convolution    of    up    to    four    exponential 
random   variables   in   the   program.    The    formulas   that   were   used 
are    given    in    Appendix    A. 
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Two  programs  are  provided  in  this  paper. 

Program  1  can  be  used  when  all  rates  are  dissimilar  cr 
all  are  the  same.  It  uses  the  formulas  on  pages  37  and  38 
only. 

Program  2  can  be  used  fcr  the  general  case.  It  makes  use 
of  all  the  formulas  givan  in  Appendix  A.  The  program 
includes  a  sorting  routine  that  determines  the  applicable 
formula  from  the  entered  failure  rates. 

A  user  guide  to  the  two  programs  is  provided  in  Appendix 
B. 
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V.        SUMMARY 
Thare   is   a    reliability  shorthand   that    denotes   the    survival 
function   of  a    system,    assuming   that    the    failure   rates    of   all 
components   are   constant. 

This   shorthand    can   be   implemented    on    the   TI-59    handheld 
calculator.    With   failure   rates,    time   to    failure   and   short- 
hand   as   input   the    TI-59    calculates   the    survival    probability 
of   the   system. 

Knowledge   of   calculus   is   not   nac?ssary  to    use    this 
method,    whereas    the   standard   procedure,    finding   the    survival 
probability  by   convolution,    requires    knowledge   of    integral 
calculus. 

The    choice    of   the    TI-59    as    the   computing  .device    for   the 
implementation    of    the    shorthand,    though,    impliad    limita- 
tions;   the    number    of   failure   rates   is    limited   due    to   the 
limited    storage    capacity    of    tha    TI-59,    and    computing    times 
are    comparatively    long.    Tha    TI-59    can    therefore   only    be   used 
for    smaller   systems,    prafarably    for   tha    solution    of   class- 
room   problems. 

For    the   solution   of   larger    problems,    the   shorthand 
should   be   implemented    on    a    state-of-the-art    personal 
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computer   using   a   general    algorithm    for    the   convolution    of 
any    number    of  exponential   random    variables. 
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APPENDIX    A 


CONVOLUTION    FORMULAS 


Appendix   A   contains    formulas   for    the    convolution    of   up 
to    four   exponential   random    variables. 

For   the  two    special   cases,    when    ill    random    variables 
have   the  same    failure   rata    and   all   have   different    failure 
rates,    general   formulas  for    the    convolution   of  any   number   of 
exponential  random   variables  are   giv  =  n. 

.These   formulas    are   ussd    in    the    two    TI-59   programs    pro- 
vided   in   Appendix   B. 
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System: 


Shorthand:  EXP(/i) 


Survival   Function:    F(t)    =   e 


-X,i 


Description : 

A    single   active   ccmponsnt   with   csnstant    failure    rate   A, 
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System: 


*\ 


s  s 


Shorthand:  EXP(  A   )     ♦    EXP  (  I   )     +    ...    +    EXP  (  A   ) 


(I*)0  (U)' 


(A*)  -It 


Survival   Function:    F(t)    = 


01 


A\ 


+  •  ♦ '  + 


(hHl! 


e 


(li)U~A     -At 


ii-AV. 


I  ~4 


Description : 

A    single   active   component   with   constant    failure    rate    is 
supported   by   n-1    identical    spares. 
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System: 


Shorthand: 


EXV{iA)     +    EXP(>U)     +    ...     +    EXP(Xu) 


Survival   Function:    F(t)    = 


^i  ~^l^ 


Is    -  A>i 


i=4        i+i 


Description: 

A   single   active   component   with   constant    failure   rate    is 
supported    by   n- 1    spares.    The   active    component   and    the    spares 
have      all   constant,    but    dissimilar    failure    rates. 
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System: 


S,  sa 


Shorthand: 


EXP(^)      +    EX?U2  )     ♦    EXP(i2  ) 


Survival  Function:    F(t)    =   Ae  +    (B    +    Ct)    e 

where      A   = 


iz 


B    =    1     -    A 

A/a  Az 


C    = 


Description: 

A  single  active  component  with  constant  failure  rate  /i^ 
is  supported  by  two  spares  with  identical  constant  failure 
rate  /i2  . 
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System: 


Shorthand: 


EXP(0      +    EXP(/i2)     +    EXP  (A,)     +    EXP(/L) 


Survival   Function:    F(t)    =   Ae     "       +  (B    +    Ct   +   Dt   )  e 

where      A  = ; 

B  =     1     -    A 

C  =    iz   - 


X; 


U„     "   lz) 


D    = 


^i/.  L 


A      ^Z 


I  {  A*  -Az) 


Description : 

A  single  active  component  with  constant  failure  rate  A,, 
is  supported  by  three  spaces  with  identical  constant  failure 
rate  A2 . 
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System: 


^ 


Shorthand: 


EXP(/l„)      +    EXP(^)      ♦    EXP(/tj)     +    EXP  Uz  ) 


Survival   Function:    P(t)    =    (A   +    Bt)  e  +     (C   +    D*)e 


where      A   = 


B   = 


A*z      ""*     «5  /v z       "** 

A**     A>z 


C  At-/*)* 


C    =    1    -    A 


D    = 


A*      A,z 

U,-lz) 


Description: 

A  single  active  component  with  constant  failure  rata  A,* 
is  supported  by  one  identical  spare  and  two  spares  with  dis- 
similar, constant  failure  rate  /t2  . 
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System: 


Shorthand: 


EXP(yi„)      +    EXP(/{2)     +    EXP(/i5)     +    EXP(/tj) 


Survival   Function:    F(t)    =   Ae  ♦   B=     z        +    (C   +    Dt)  e      5 


where      A    = 


A  z  /ti 


z 


B     = 


C     = 


D    = 


A.4     Ar% 


a*  At ,  A*kz  k^  (     'i 


A 


A*    kz   A.3 

(  A,«-A>i)  (4,z  "^3) 


Description: 

A   single    active   component    with   constant    failure   rate     A*» 
has    three   spares.    One    spara    has    constant    failurs    rate    <iz  , 
two    spares    are    identical   with  constant    failure   rate    A.3  . 
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APPENDIX    B 
USEE    GUIDE    TO    TI-59    PROGRAMS 

Appendix   B    contains  a   user   guide    to   two   TI-59    programs, 
which    use    reliability    shorthand    and   failure   rates    as    input 
to   compute   the   survival    probability   of    a    system. 

PROGRAM    1    is    designed   for   the   two    special   cases   where 
the    reliability   shorthand   is  of  the   form 

EXP  (  K    )     *    EXP  (  X    )     +    ...     f    EXP (  X    ) 
or 

EXP(/i„)     +    EXP(i2)     +    ...     f    EXP  (  /£,,). 
In   the   first    case   the    number   of   terms    is   not    limited,    wher- 
eas   in   the    second    case  the    number   of    terms    is    limited    to    40 
due    to   limited   storage  capacity   of   the    TI-59.    In    this   case 
the    number    of   terms   can    be    increased    to    70    by   entering    9    in 
the    display  and   pressing    2nd   Op    1    7    . 

PROGRAM    2    is   designed   to   solve   problems   of   the    kind, 
that    were    introduced    in    Chapter    III. 5.     .    Due    to    limited 
memory   of    the   TI-59    the   number   of   exponential    terms   under 
one    weight    in    shorthand   notation    is    limited   to   four. 
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All   results    will   be    printed,    if   the    TI-59   is    connected 
to    a    TI    FC-100A   or   TI    PC-10  0C  printer. 
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PROGRAM    1    :    Procedure 

1.  Use   any    library    module.     Read   in   program    1     (side    1    of 
the   magnetic   card) 

2.  Press    2nd    C    to    initialize. 

3.  Enter   n,    the    number  of   exponential    terms   to   be    convolved, 
in   the   display  and    press    A. 

4.  Enter  time   t    and  press   B. 

5.  Enter  h, [  and    press   C    .    When   all    failure    rates    are   -he 
same,    enter  /i   only   once. 

6.  a)    To    find   the    survival    probability   of    the   system, 

when    all    failure  rates   are  the   same,    press    2nd    A'. 
b)    To   find   the    survival    probability    of    the   system, 

when    all    failure  rates  are   dissimilar,    press   2nd   B" . 
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PROGRAM  1  :  Sample  Problems 

1.  Find  the  survival  probability  of  a  parallel  system 

(  compare  Chapter  III. A. 2  ) 

a)/i  =  .3   ,   t  =  7   ,   n  =  2 

b)    Shorthand   notation: 
EXP(.6)     +    EXP(.3) 


c) 


Enter     Comment    Prass     Display 
Initialize     C         0 


2 

n 

7 

t 

.6 

2l 

.3 

I 

\ 


5 


F(t) 


7 

.3 

.  3 

.2299172797 


calculation  takes  13  seconds 
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2.    Find   the  survival    probability  of   a    standby-system   with 
dissimilar   components    (    compare   Chapter   III. A. 3   )     . 

a)  iA  =    .4      ,  /i2=    .5      ,      t    =   6  ,      n    =    2 

b)  Shorthand    notation: 

EXP(.4)     +    EXP(.5) 

c)  Enter             Comment  Prsss             Display 

Initialize  Z%                      0 

2                          n  A                          0 

6                           t  3                           6 

.4                        X4  C                           .4 

.5                        /i2  C                           .5 

F(t)  B'                          .254441493 
calculation    takes    13    seconds 
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3.    Find  the  survival    probability  of   a    standby-system    with 
one   active   component   and    four   similar   spares. 

a)  /(,  =    .  3      ,      t=7      ,      n   =   5 

b)  Shorthand    notation: 

EX?(.3)     ♦    EXP(.3)      ♦    EXP(.3)     +    EXP(.3)     +    EXP(.3) 

c)  Enter     Comment   Pr?ss     Display 

Initialize     Z%  0 

5  n        A  0 

7  t         B  7 

.3  i  Z  .3 

F(t)      A*  .93787388U8 

calculation  -akes  9  seconds 
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PROGRAM    2    :    Procedure 

CASE    I    :    To   find   the    convolution    of    up    to    four    exponential 
random   variables. 

1.  Use   any    library    module. 

Re-Partition    (    enter    2   in   the    display,    press   2nd    Op    17    )  . 
Read    in    all    four    sides  of   the    magnetic    card. 

2.  Press    2nd   C    to    initialize. 

3.  Enter   n,    the    number   of   exponential    terms   to   be    convolved, 
in   the   display  and    press    A. 

4.  Enter  time  t   and   press   B. 

5.  Enter  L^    and   press    C    (   n    entries   )     . 

REMARK:    Failure    rates,    which    appear    only   once    in    the 
expression,    haze    to    be    entered    before    failure 
rates,    that  appear   several    times. 

6.  To    find    the   survival   probability   of    the    system    press   E. 
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PROGRAM    2,    CASE    I    :    Sampls    Problems 

(1)     Shorthand   notation 

EXP(/i„)     +    EXP{/i2)     +    EXP(/te) 

Sample    values   :  /i4    -   .3    ,   A^   -    .4    ,      t   =    7 

Procedure    : 

Enter     Comment    Press      Display 
Initialize     Z%  0 

0 
7 

.3 
.4 
.4 
F(t)  E  .5363473366 

calculation    takes    14    seconds 


3 

n 

i 

7 

t 

3 

.3 

L, 

C 

.4 

iz 

w 

.4 

li 
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(2)     Shorthand   notation 

EXP(yt„  )     «■    EX?{il  )     +    EXP(/t2)  +    EXP  (  X2  ) 

Sample   values   :   /£,,    =  .2    ,    A>i  =  ,4    ,      t   =    3 
Procedure    : 

Enter             Comment        Prsss  Display 

Initialize             2*  0 

A  0 

B  3 

Z  .2 

C  .4 

c  .a 

c  .4 

E  .9809746099 
calculation    takes    20    seconds 


4 

n 

3 

t 

.2 

Kj\ 

.4 

It 

.4 

li 

.4 

iz 

P(t) 
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(3)     Shorthand  notation 

EXP  ( yt„   )     +    EXP  (  lA  )     +    EXP  (  it  )  +    EXP  (  A,t ) 

Sample   values   :   A>+    -   .4    ,    Az  =  .  3    ,      t    =    5 
Procedure    : 

Enter             Comment        Pr  =  ss  Display 

Initialize             C '  0 

A.  0 

B  5 

C  .4 

C  .3 

C  .3 

E  .9029040721 
calculation  rakes  20  seconds 


a 

n 

5 

t 

.4 

I* 

.4 

Ka 

.3 

it 

.3 

it 

?(t) 
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(4)     Shorthand   notation 

EXP  (/{,„)     +    EXP(/£2)     +    EXP{/3)  +    EXP(/i3) 

Sample    values    :   /i^    =    .1    ,     hz  =  .3    f    /lj  =    .5    , 

t    =    10 

Procedure    : 

Enter             Comment        Pr=ss  Display 

Initialize             3'  3 

A  0 

B  10 

C  .  1 

Z  .3 

C  .5 

C  .5 

F  (t)               E  .7312684703 
calculation    takes    25    seconds 


4 

n 

10 

t 

.1 

I, 

.3 

hi 

.5 

/1>3 

.5 

ii 
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PROGRAM    2    :    Procedure 

CASE   II    :    to   solve   problems    of   the   kind,    that    were 
introduced    in   Chapter   III.  3.    . 

1.  Derive   the   system's   shorthand   notation.    Find   either    the 

-  graphical   representation   or 

-  the    MIX-notation    . 

2.  Use   any    library    module. 

Re-Partition    (   enter    2   in   the    display,    press   2nd   Op    17   ). 
Read    in    all    four   sides   of   the    magnetic    card. 

3.  Press   2nd   C    to   initialize. 

4.  Enter   time   t    and  press   3. 

5.  Repeat   the    following   steps   for    each    path    of   the    graphical 
representation,    i.e.    for    each   convolution   in   the 
MIX-ncta-ion. 

a)  Enter   n,    the    number   of  exponential   terms   to    be 
convolved,    in   the    display    and    press    A. 

b)  Enter    X^    and   press   C. 

REMARK:    Failure    rat.es,    which    appear   only   once    in   the 
expression,    have   to    be    entered    before   failure 
rates,    that    appear   several   times. 

c)  Enter   p« f    the   weight    in  the   ith    path,    and   press    D. 

d)  Tc    find   the    part   of   the  system's    survival    probability, 
that    is   contributed   by   the    ith    path,    press    E. 
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6.    To    find    the    survival    probability   Df   the   system 
press   2nd    E • . 
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PROGRAM    2,    CASE   II    :    Sampla    Problems 

1.    Find   the   survival    probability   of   a    parallel  system 

with    dissimilar    failure    rates    (compare    Chapter    III.B.1) 

a)  LA  -    -1      ,     J>z    =   .2      ,      t    =   2 

b)  Shorthand   notation 


EXP  (.  1) 


(.2/.3) 


EXP(.3)     + 


(-1/.3) 


EXP  (.  2) 


F(t)     =    MIX[      (.2/. 3)  (     EXP(.1)     ♦    EXP  (.3), 
(.  1/.3)  (     EXP(.  2)     +    EXP  (.3)   ] 
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c) 


Procedure  : 

Enter     Comment  Prsss  Display 

Initialize  Z*                       0 

2          t  B          2 


2 

ni 

A 

0 

.1 

*« 

Z 

.1 

.3 

/I,,  +  Ki 

3 

.3 

(.2/. 3) 

P, 

D 

.6666666667 

E 

.635793541 

2 

nE 

A 

0 

.2 

4» 

.2 

.3 

/V/J     +      /Cl 

r^ 

.3 

(.  V.3) 

p2 

D 

.3333333333 

E 

.304445622 

F(t) 

E1 

.940239163 
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2.  Find  the  survival  probability  of  a  series  system  with 
one  spare  as  introduced  in  Chapter  III.B.2  . 

a)  /m  =  .3   ,  A*2   =  -5   ,   t  =  7 

b)  Shorthand    notation 

ZEPO 
(-3/.8) 

EXP(.8)     + 


(-5/.8) 


EXP  (.  8) 


F(t)     =    MIX[      (.3/. 8)  (     EXP(.8), 

(.5/. 8)  (    EXP(.8)     +    EXP(.8)], 
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c) 


Procedure    : 


Enter 

Com  mer.t 

Prsss 

Display 

Initialize 

C  ' 

0 

7 

4- 

B 

7 

1 

n„ 

h 

0 

.8 

'Vl    +   Kx 

n 

.3 

(.3/. 8) 

p-. 

D 

.375 

E 

.0013866989 

2 

n2 

A 

0 

.8 

/t-i  +   /t^ 

3 

.8 

.8 

/C-,  +    ^z. 

.8 

(.5/. 8) 

?z 

D 

.625 

? 

.0  152536878 

F(t) 

S' 

.0166U03867 
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3.    Find   the   survival   probability   of   a    Two-out-of-Three 
System   as    introduced   in    Chapter    III.B.3    . 

a)  KA  =    -2      ,      iz    -  .U      ,      /i3  =    .5       ,      t    =    9 

b)  Shorthand   notation 

EXP  (.  9) 
(.2/1.1) 


(.4/1.1) 


EXP(1.1)     + 


EX?  (.  7) 


(.5/1.1) 


EXP  (.6) 


F(t)  =  MIX[  (.2/1.1)  (  EXP(.9)  +  EXP(1.1)), 
(.4/1.1)  (  SXP(.7)  +  EXP(1.1)), 
(.5/1.1)   (    EXP(.6)     +    SXP(1.1))] 
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c) 


Procedure    : 


Enter 


Comment        Press 


Initialize 


(.4/1.1) 


2 
1.1 
.6 
(.5/1.1) 


A>*    +    4*2, 


r>  t 


2                          n„ 

A 

1.1               A*  +     A*2   +  A,$ 

r* 

.9                           ^2    +    /ij 

r* 

(-2/1.1)                 p^ 

D 

E 

2  n2  A 

1.1         /l^  +   Ai  +  A3       c 

.7  /(,,    +     /I3  c 


F(t) 


D 

E 


Display 

0 

9 

0 

1.1 

.9 

. 1818181818 
.0002624871 

0 

1.1 

.7 

.3636363636 
.0018043754 

0 

1.1 

.6 

.U545454545 
.0044892129 
.0065560755 
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COMPUTES    LISTINGS 


PROGRAM     1 
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PROGRAM    1    continued 


L^BEL    ADRESSES 


h  Si 
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PROGRAM    2 


crp 

O   =   Li 
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PROGRAM    2    continued 


bU 


A       CM  M 

q        1  Q 
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PROGRAM    2   continued 
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PROGRAM    2    continued 


67 


PROGRAM    2   continued 
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PROGRAM    2    continued 


fiff4 
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PROGRAM  2  continued 


LABEL  ADRSSSES 


kj! 
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